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ABSTRACT 



(N ■ 

■ Aims. We aim at estimating the spectral energy distributions (SEDs) and the physical parameters related to the black holes harbored in eight 
£^ ' high X-ray-to-optical (F x /Fr >10) obscured quasars at z > 0.9 selected in the 2-10 keV band from the HELLAS2XMM survey. 

Methods. We use IRAC and MIPS 24 fim observations, along with optical and A'j-band photometry, to obtain the SEDs of the sources. The 

■ observed SEDs are modeled using a combination of an elliptical template and torus emission (using the phenomenological templates of Silva 
et al. 2004) for six sources associated with passive galaxies; for two point-like sources, the empirical SEDs of red quasars are adopted. The 
bolometric luminosities and the M BH - L K relation are used to provide an estimate of the masses and Eddington ratios of the black holes 
residing in these AGN. 

\ Results. All of our sources are detected in the IRAC and MIPS (at 24 yum) bands. The SED modeling described above is in good agreement 
• • . with the observed near- and mid-infrared data. The derived bolometric luminosities are in the range « 10 45 - 10 47 erg s , and the median 
' 2-10 keV bolometric correction is w25, consistent with the widely adopted value derived by Elvis et al. (1994). For the objects with 
k> elliptical-like profiles in the K s band, we derive high stellar masses (0.8 - 6.2)xl0" M Q , black hole masses in the range (0.2 - 2.5)xl0 9 M Q , 
t , ■ and Eddington ratios L/Lem < 0.1, suggesting a low-accretion phase. 
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1. Introduction 

Hard X-ray surveys have clearly revealed the important role 
played by obscured Active Galactic Nuclei ( AGN) to repro 



duce the cosmic X-ray background (XRB; e.g. JComastri et al 



1995) and have provided evidence that a significant fraction of 
the accretion-driven energy density in the Universe resides in 



obscured X-ray sources (e.g., Bars ler et alj2005l Hopkins et al 
2006:lHickox & Markevitch 200i 



Until recently, the limited information on the broad-band 
emission of the counterparts of obscured X-ray sources pre- 
vented a reliable determination of their bolometric luminosi- 
ties. The lack of a proper knowledge of the spectral energy dis- 
tributions (SEDs) of obscured sources has led many authors to 
adopt, in the computatio n of the bol o metr ic luminosities, the 
average value derived by Elvis et alj d 19941) . although in that 



work the sample comprises mostly local unabsorbed quasars. 
By the current work, we aim at providing a robust estimate of 
the bolometric luminosity for obscured sources, which is an 
essential parameter to derive the cosmic mass den sity of super- 
massive black holes (SMBHs, i.e., following the Soltan 1982 
approach). A reliable estimate of the bolometric luminosity of 
obscured AGN is typically limited by the actual capabilities 
of disentangling the nuclear emission (related to the accretion 
processes) from that of the host galaxy which, unlike for un- 
absorbed quasars, often dominates at optical and near-infrared 
(near-IR) bands. 

A significant fraction of high-redshift, luminous obscured 
AGN (the so-called Type 2 quasars) may have escaped spec- 
troscopic identification due to their faint optical counterparts, 
thus preventing current studies from an accurate sampling of 
obscured sources. Mid-infrared (mid-IR) observations appear 
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Table 1. Properties of our targets 



Source Id. 


2- 

(io- 


-10 keV flux 1 " 

14 -2 - 

erg cm s 


R 

') 


K s 


Morphea) 


X/O 


z 


Abell 2690#75 




3.30 


24.60 


18.33 


E 


1.86 


2.13" 


PKS 0312-77#36 




1.90 


24.70 


19.13 


E 


1.66 




PKS 0537-28*91 




4.20 


23.70 


18.99 


E 


1.60 




PKS 0537-28*54 




2.10 


25.10 


18.91 


E 


1.86 




PKS 0537-28*111 




2.10 


24.50 


17.64 


E 


1.62 




Abell 2690*29 




2.80 


25.10 


17.67 


P 


1.99 


2.08 i 


PKS 0312-77*45 




2.80 


24.40 


18.62 


P 


1.71 




BPM 16274*69 




2.27 


24.08 


17.87 


E 


1.49 


1.35* 



t 2-10 kev X-ray fluxes from lPerola etai] d2004l) . 

" Tentative spectroscopic redshift from near-IR spectroscopic obs ervations iMaiolino et al]2 006). 



Spectroscopic redshift from near-IR spectroscopic observations (Ma iolino et al.i l20061 



to be fundamental for this class of objects, since they are only 
marginally affected by dust obscuration and are able to recover 
the nuclear emission. With mid-IR observations, we expect to 
reveal the nuclear radiation re-processed by the torus of the 
obscured active nuclei, which are often recognized as such by 
means of their X-ray emission only. For these sources, the 
soft X-ray emission, which is photo-electrically absorbed by 
the gas, and the optical emission, extincted by the dusty cir- 
cumnuclear medium, are expected to be downgraded in energy 
and to emerge as thermally reprocessed radiation in the IR at 
wavelengths in the r ange between a few and a few hundred /im 
dGranatoet al.lll997l) . 

The potentialities of m id-IR observations were firstly 
shown by iFadda et all d2002l) . who detected with ISOCAM at 
15 yum about two-thirds of the X-ray sources detected in the 
5-10 keV band in the XMM-Newton Lockman Hole survey. A 
simil ar high de t ection rate at 24 im\ has been recently reported 
by iRigbv et alj d2004l) and iFranceschini et alj ([2005) studying 
the Spitzer counterparts of the Chandra sources in the CDF- 
S and i n the ELAIS-N1 field , respectively, within the SWIRE 
survey ( Lonsdale et alJl2004 ). 

In this context, a new interesting class of objects is emerg- 
ing from the current X-ray surveys: these sources are character- 
ized by a high (> 1 ) X-ray-to-optical flux ratio (hereafter X/ojH 
for comparison, unobscured Type 1 AGN have a broad distri- 
bution peaked at X/O-0. Objects with X/O >1 are about 20% 
of the hard X-ray selected sources, and the fraction of these 
sour ces seems to remain con stant over « 3 decades of X-ray 
flux ( Comastri & Fiorell2004 ). By definition, sources with high 
X/O are among the faintest sources in the optical band. In the 
shallow, large-area X-ray surveys (e.g., the HELLAS2XMM 



survey, with F2- 10 kev > 10 erg cm s over « 3 deg ; 
Baldi et al. 2002), where the identification of a sizable sam- 



ple of sources with X/0>1 has been possible (e.g.. lFiore et al.1 
120031) . the X/O selection criterion has proven to be effective in 
selecting Type 2 quasars at high redshifts. 



We have performed a pilot program to study with Spitzer a 
sample of eight sources selected in the 2-10 keV band from the 
HELLAS2XMM survey on the basis of their high (> 1) X/O and 
large column densities (Nu >10 22 cirT 2 ). The sample observed 
with Spitzer has been previously investigated in other bands. 
The most surprising finding of the follow-up campaigns was 
the association of th ese sources with luminous near-IR objects 
(IMignoli et al ] |2004l) . placing them into the class of Extremely 
Red Objects (EROs, R - K >5). 

The outline of the paper is as follows: in Sect. 2 we present 
the sample selection; Spitzer data reduction and analysis are 
discussed in Sect. 3, while in Sect. 4 we describe the analy- 
sis of the SEDs. Finally, in Sect. 5 we estimate the bolometric 
luminosities, the stellar masses of the host galaxies, the black 
hole masses, and the Eddington ratios. 

Throughout this paper we adopt the "concordance" 
(WMAP) cos mology (Hg=7 km s 1 Mpc~\ O M =0.3, and 
Qa =0.7: fepergel et aLll2003l) . Magnitudes are expressed in the 
Vega system. 



2. Sample selection 

The eight objects presented in this paper (see Table [TJ were se- 
lected among the 10 HELLAS2XMM high X/O ratio sources 
detected in the K s band with ISAAC at ESO-VLT; for details 
on the asso ciation of the A^-ba nd counterpart to the X-ray 

d2004T) . Two sources of the origi- 



Migno li et al 



source, see 

nal sample were not selected for Spitzer observations: one 
(PKS 0537-28#31) is associated with a disky galaxy, while 
the other (BPM 16274#181) has an ambiguous A^-band mor- 
phological classification. All but one of the sources observed 
by Spitzer belong to the firs t square degree field (122 X-ray 
sources; see lFiore et alj|2003 for the spectr oscopic and photo 



metric identification and Perola etal]|2004 for the X-ray spec- 
tral analysis); the only exception is BPM 16274#69, which be- 
longs to the sample of the secon d square degree (110 X-ray 



1 X/O is defined as log 



log F x + £+5.5. 



sources; see 
tion). 



Cocchia et al.l 12007; Lanzuisi et al., in prepara- 
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Abell 2690#75 PKS 0312-77#36 PKS 0537-28^91 




PKS 0537-28^54 PKS 0537-28#1 1 1 Abell 2690#29 




PKS 0312-77#45 BPM 16274#69 




Fig. 1. ISAAC K s images, centered on the K s counterpart of the X-ray sources; each box is 30" wide. North is to the top and East 
to the left. Contour levels of the 24 yum emission corresponding to [3, 4, 5, 7, 10, 20, 40]cr are superimposed to each image. 



The selected sources, although faint in the optical band 
(23.7<R<25.1), are all bright in the near-IR ( K<<\9) and, all 
but on e, have R - K s > 5, thus being EROs. Mi gnoli et al.1 
(2004) were able to study the surface brightness profiles of 
these sources in the K s band and obtained a morphological clas- 
sification. Only two sources are classified as point-like objects, 
while all of the others are extended, with clear detection of the 
host galaxy and radial profiles consistent with those of ellipti- 
cal galaxies. In this latter class of sources, the central AGN is 
evident in the X-ray band, while in the near-IR the host galaxy 
dominates. An upper limit to the contribution of a central un- 
resolved source (i.e., the nuclear emission), ranging from 2% 
to 12% of the galaxy emission, was obtained. Furthermore, us- 
ing both the R - K colour and the morphological information, 
a minimum photometric redshift, in the range 0.9-2.4, was es- 
timated for these sources. 



For three of the sources with the reddest colours 
(Abell 2690#75, Abell 2690#29 and BPM 16274#69; see 
Table Q}, near-IR spectroscop ic observations with I SAAC at 
ESO- VLT were performed by iMaiolino et al. (2006), thus al- 



lowing for a spectroscopic identification for at least two of 
these sources (one redshift measurement appears tentative). 
The point-like source (Abell 2690#29) shows the typical rest- 
frame optical spectrum of high-redshift du st-reddened quasars, 
with a broad Ha line dGreggetal]l2002h . The other two ob- 
served sources, both of them extended in the K s band, have nar- 
row emission-line spectra: one is a LINER-like object at z= 1 -35 
and the second source has a spectrum with a single weak line, 
tentatively associated with Ha at z=2.13. Consistently with the 
morphological information, in the first source the AGN domi- 
nates the emission, while in the other two sources the nuclear 
spectrum is heavily diluted by the host galaxy starlight. 
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Table 2. Spitzer flux densities 



Source Id. 


3.6 pm 


4.5 yum 


5.8 yum 


8.0 yum 


24 yum 




S V ±AS„ 


S V ±AS„ 


S V ±AS„ 


S V ±AS V 


S V ±AS V 


Abell 2690#75 


51 ±5 


56 ±6 


89 ±11 


139 ± 15 


565 ± 62 


PKS 0312-77*36 


41+4 


44 ±5 


40 ±8 


71+9 


236 ± 30° 


PKS 0537-28*91 


28+4 


35 ±4 


42 ± 8 


80 ± 10 


301 ± 40 


PKS 0537-28*54 


31 ±4 


35 ±4 


50 ± 10 


47 ± 8 


279 + 45 


PKS 0537-28*111 


88 ±9 


75 ± 8 


41 ±6 


46 ±7 


148 ± 28 


Abell 2690*29 


141 ± 14 


185 ± 19 


260 ± 27 


371 ±38 


1012 ±106" 


PKS 0312-77*45 


50 ±6 


62 ± 7 


69 ± 10 


78 ± 10 


249 ± 35 


BPM 16274*69 


86 ±9 


92 ±9 


97 ± 11 


120 ± 13 


286 ± 34 



The flux density is reported in units of yijy. " The 24 yum flux density is probably over- 
estimated due to contamination from nearby sources and should be considered as an 
upper limit. 



3. Spitzer observations and data reduction 

The whole sample of eight hard X-ray selecte d sources has 



been observed by Spitzer (IWerner et al.l 120041) . with IRAC 



(Fazio et al. 2004) observations of 480 s integration time and 
MIPS dRiekeetal.ll2004l) observations at 24 yum with a total 
integration time of « 1400 s per position. IRAC observations 
were performed in photometry mode with frame time of 30 s 
and dither pattern of 16 points. The MIPS 24 pm observations 
were performed in MIPS photometry mode with frame time of 
10 s, 10 cycles and small-field pattern. To reduce overheads, 
the cluster option was used when possible. 

For the IRAC bands, we used the final combined post-basic 
calibrated data (BCD) mosaics produced by the Spitzer Science 
Center (SSC) pipeline (Version S12.0-S13.01). At 24 pm, 
we started the analysis from the BCD produced by the SSC 
pipeline (Version S12.4.2-S13.01) and then we applied ad hoc 
procedures to optimize the reduction, since some of our sources 
were close to the detection limit (see Table [2]). We remind that 
BCD are individual frames already corrected for dark, flat field 
and geometric distortion. We improved the quality of the BCD 
by correcting each individual BCD fo r a residual flat fie ld de- 
pending on the scan mirror position dFadda et al. 2006). The 
residual flat fielding was obtained from our own data by aver- 
aging the BCD corresponding to the same scan mirror position 
and the same Astronomical Observation Request (AOR), con- 
sidering all the different cluster positions. To each BCD, its 
median level was subtracted before this operation. This pro- 
cedure was possible since our observations are not dominated 
by background fluctuations. The corrected BCD were co-added 
and background-subtracted using the SSC MOPEX software 
(Ma kovoz & Marleaul2005l) . The resulting mosaics were made 
with 2.4" pixel size. The overall analysis at 24 pm produces 
mosaics with a typical noise of « 0.020 MJy/pixel (a factor of 
2 lower in comparison with the SSC pipeline mosaics). The 
noise map has been computed for each mosaic by scaling the 
measured mean rms of the central part of the map according to 
the inverse square root of the coverage map. 

The flux densities of our targets in IRAC and MIPS bands 
were measured on the signal maps using aperture photometry at 
the position of the sources. The chosen aperture radius for the 



IRAC bands is 2.45" and the adopted factors for the aperture 
corrections are 1.21, 1.23, 1.38 and 1.58 (following the IRAC 
Data Handbook, Version 3.0) at 3.6 pm, 4.5 yum, 5.8 pm and 8 
pm, respectively. 

The chosen aperture radius at 24 pm is 7.5"; aperture correc- 
tions were derived by examining the photometry of bright stars. 
Taking into account an additional correction of 1.15 to match 
the procedure used by the MIPS instrument team to derive cal- 
ibration factors from standard star observations, the resulting 
ape rture correction is 1 .57 (in agreement with the SWIRE team, 



Sura ce et al.ll2005h . 



Table |2] reports the results of the Spitzer observations. To 
compute the photometric uncertainties, we added in quadra- 
ture the noise map and the systematic uncertainties («10%, see 
MIPS and IRAC Data Handbook 2006, Version 3.0). The rel- 
ative photometric uncertainties range from «10% in the best 
cases (IRAC channels 1 and 2), up to «20% at 24 pm in the 
worst ones. 

All of the eight sources are clearly detected in both IRAC 
and MIPS 24 /im bands. At 24 pm, the flux densities span an 
order of magnitude, ranging between «1000 yt/Jy and a: 150 pJy, 
with the faintest source (PKS 0537-28#lll) close to the 5cr 
detection level. 

In Fig. Q] the /T,-band images, along with the contour lev- 
els of the 24 pm emission, are shown. At 24 //m, the sources 
PKS 0312-77#36 and Abell 2690#29 appear to be confused. 
In particular, in both cases, there is a second source at ss 8-10", 
unrelated to the targets. The contribution of these sources to 
the 24 pm flux density of our targets has been estimated by a 
decomposition analysis (using the PSF fitting algorithm IMFIT 
within the AIPS environment). Furthermore, both the sources 
PKS 0312-77#36 and Abell 2690#29 present a second object 
at « 2", clearly visible in the K s images, too close to our targets 
for a decomposition analysis, given the 24 yum pixel size. Since 
these close companions become increasingly fainter moving 
from the K s bands to the longer IRAC wavelengths, we have 
attributed the entire flux density estimated from the decompo- 
sition analysis to our targets (see Table©. However, the 24 yizm 
flux densities for these two sources should be treated as upper 
limits (see Fig. [2]). 
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We point out that the deblending procedure measures the 
peak flux density. To convert the peak flux density into total 
flux density, we have assumed that the deblended objects are 
point-like sources and applied a correction factor of 8.9 derived 
from the 24 pm Spitzer PSF and including the 1.15 calibration 
factor. 

4. Analysis of the spectral energy distributions 

In this section we provide an analysis of the SEDs of our 
sources, in order to derive the energy distribution of the nu- 
clear component "cleaned" by the host galaxy contribution. As 
anticipated in Sect. [T| the determination of the nuclear SEDs 
over a large wavelength range is an essential step to estimate 
the physical properties of the black hole, such as its bolometric 
luminosity, mass, and accretion rate. Taking advantage of the 
new Spitzer photometric points, simultaneously to the SED de- 
termination we have estimated the photometric redshifts of our 
sources. These new values a re then compared wi th the min- 
imum redshifts estimated by IMignoli et al.1 d2004l) using only 
the R and K s ba nds and with the spect roscopic ones measured 
in three cases by lMaiolino et al.l (120061 see 

Given the different morphological properties of our 
sources, two approaches have been adopted, one for the sources 
dominated in the K s band by the host galaxy (elliptical-like 
sources) and another one for the sources dominated by the nu- 
clear component (point-like sources). 

4.1. Elliptical-like sources 

From the /f s -band morphological analysis, we know that at 
least up to the observed 2.2 pm band the stellar contribu- 
tion dominates the emission in these sources. At longer wave- 
lengths, the nuclear component is expected to arise as repro- 
cessed radiation of the primary em ission, while the stellar 



component is ex pected to drop (e.g.. lBruzual & Charlolll2003 
Silva et alJEooH) . 



In the analysis presented in this work, we have decided 
to follow a phenomenological approach, checking whether the 
emission of our sources can be reproduced as the sum of two 
components, one from the host galaxy and the other related to 
the reprocessing of the nuclear e mission by the d usty torus en- 
visaged by unification schemes (lAntonuccil 1 1 993b . The shape 
and relative strengths of the two components have to be con- 
sistent with all our observed data sets (multi-band photometry, 
/fj-band morphology and magnitude, /fj-band upper limit on 
the nuclear component, and X-ray spectral analysis). 

For the galaxy component, we adopted a set of six galaxy 
templates, obtained from the synthetic spectra of GISSEL 2003 
( Bruzual & Charlotl2003 ) assuming a simple stellar population 
and spanning a wide range of ages, from 1 Gyr up to a "maxi- 
mum age" model (z form = 20). The "m aximum age" model has 
been adopted by IMignoli et al. (2004) to derive the minimum 
photometric redshift for the sources of the current sample. 

For the nuclear co mponen t, we adopted the nuclear tem- 
plates from Silva et al. (2004J ), base d on the radiative transfer 
models of lGranato & Danese ( 1994). We chose these templates 
since in the work of ISilva et al.l (120041) the radiative transfer 



models are used to interpolate the observed nuclear IR data for 
a sizable sample of local AGN. We must note, however, that 
the nuclear observed SEDs are available only in the 2-20 pm 
regime, where data from small-aperture instruments are avail- 
able. At wavelengths above « 20 pm, the SEDs are model ex- 
trapolations^ 



Silva et al.l (120041) found that the nuclear SEDs can be ex- 



pressed as a function of two parameters, the hard X-ray (2- 
10 keV) intrinsic luminosity, which provides the normalization 
to the SED, and the colu mn density Nh, wh ich gives the shape 
to the SED (see Fig. 2 in lSilvaet alll2004 . The shapes of the 
SEDs of the Seyfert galaxies are assumed to be valid also at 
quasar luminosities. 

In the attempt to provide a better estimate for the source 
redshifts, we us ed the four torus templates as given by 



Silva et al.l (120041) . which depend on the column densities Nh, 



and we left the normalizations free to vary. The redshift interval 
explored by this procedure is 0.5-3.0; the best-fit solution is ob- 
tained when the algorithm, based on the x 2 , finds a minimum 
in the galaxy template, torus template and redshift parameter 
space. 

For four out of six sources the procedure finds a clear min- 
imum x 2 ' which allows us to determine a photometric red- 
shift with relatively good accuracy (see Table |3). For sources 
BPM 16274#69 and PKS 0537-28#54, our procedure con- 
strains only the lower bound of the redshift interval. For all 
the six sources, the estim ated redshifts are consistent with 
the minimum redshifts of IMignoli et al.l (120040 . In case of the 
source BPM 16274#69, where a secure spectroscopic redshift 
is available, the minimum photometric redshift {z p hoi > 1-25) 
is consistent with the spectroscopic one (z = 1.35). For source 
Abell 2690#75, we find z pho , = 1.30+°;|g, which is signifi- 
cantl y lower than the sp ectroscopic value (z = 2.13) reported 
by Maiolino et al. d2006h . In this case, we choose to adopt the 
photometric determination, since the spectroscopic redshift is 
based on the tentative detection of a single line. 

The results of the SED fitting and decomposition are shown 
in Fig. |2] The dot-dashed line represents the best-fit galaxy tem- 
plate, the dashed line is the best-fit nuclear template and the 
thick solid line is the sum of the two components. The best-fit 
galaxy templates are all typical of early-type galaxies with ages 
between 3 and 6 Gyr. We find an overall agreement between the 
SED templates and the data points. Moreover, in all but one of 
the sources, the nuclear component derived from the best fit 
is consistent with the upper limits derived from the analysis 
of the A'j-band images (shown as downward-pointing arrows). 
In Fig. [2] we also report as dotted line the SED of the nuclear 
component normalized to the intrinsic (i.e., de-absorbed) X-ray 
luminosity following the prescriptions of ISilva et al .1 (120041) . 
where, at a given Nh, the normalization depends only on the in- 
trinsic hard X-ray luminosity. The overall agreement between 
the SEDs normalized to the X-ray luminosity and the best-fit 
SEDs is extremely interesting, being consistent within a factor 
of « 2-3. 

In Table [3] we report, along with the photometric redshifts, 
the column densities N# and the de-absorbed L2-10 keV lumi- 
nosities. We derive rest-frame N# column densities in the range 
10 22 — 10 23 ' 4 cirT 2 and 2-10 keV luminosities in the range 
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A2690#75 zP hot =1.3 



pks0312-77#36 z P hot =0.9 




1000.00 



100.00 



100.0 



[rest] 



pks0537-28#91 Z P hot =1.3 




100.0 



[rest] 



pks0537-28#1 1 1 z P hot = 1 .2 





1000.00 



100.00 



100.0 



100.0 



[rest] 



pks0537-28#54 Z P hot >1.3 




100.0 



A[/im] 



[rest] 



bpm16274#69 z s =1.35 




100.0 



[rest] 



[rest] 



Fig. 2. Rest-frame SEDs of the elliptical sources (black filled circles) compared with the best-fit model obtained as the sum (solid 
line) of an early-type galaxy (dot-dashed line) and a nuclear component (dashed line). For comparison, the nuclear component 
as derived from the X-ray normalization is also reported (dotted line) . The nuclear A", -ban d upper limits (downward-pointing 
arrows) were derived from the morphological analysis carried out by iMignoli et alJ (12004). The 24 //m upper limit of source 
PKS 0312-77#36 takes into account a possible contribution from a companion source. z s means that the redshift is spectroscopic 
(see $2]for details), while z phor means that the redshift is photometric (see 34.11 ). 



10 —10 erg s , placing these sources among the Type 2 
quasar population. 



4.2. Point-like sources 

From the /T s -band morp hological analysis, Abell 2690#29 and 
PKS 0312-77#45 (see IMignoli et alj|2004l and Table HJ show 



a completely different appearance at 2.2 //m in comparison 
with the sample of extended objects; these two sources have 
their near-IR emission mostly dominated by an unresolved 
source. The dominant role played by the AGN is supported 
for Abell 2690#29 also by its nea r-IR spectrum, where a broad 
Ha emission line is detected fsee lMaiolino et al.ll2.Q06l and $2}. 
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10000.0 




100.0 



J [rest] 



[rest] 



Fig. 3. Rest-frame SEDs of the two point-like sources (black filled circles) compared to an extinguished quasar template (dashed 
line) and th e best-fit red quasa r template (solid line). The extinguished quasar template is obtained from the unobscured quasar 
template of lElvis et al. ( 1994b using the SMC extinction law with E(B - V) = 0.7 and scaled to fit the R - K s colour (for 



comparison, the unobscured quasar template is also shown as a dotted line). The red quasar template is taken from lPolletta et al 



(120061) . The 24 pm flux density upper limit for the source Abell 2690#29 takes into account a possible contribution from a 
companion source. z s means that the redshift is spectroscopic (see §|2]for details), while z pho ' means that the redshift is photometric 
(see JO). 



Unfortunately, the near-IR spectroscopy information is absent 
for PKS 03 12-77#45 (see TableQ). 

We first tried to reproduce their observed SEDs by redden- 
in g the composite t emplate spectrum of bright Type 1 quasars 
of Elvis et al.l (119941) with several extinction laws . Reddening 
has been applied as prescri bed by ICalzettil (119971) for a dust- 
screen model and by iPeil d 19921) for the Small Magellanic 
Cloud (SMC) galaxy. The two prescriptions produce similar 
effects at A > 0.5 pm, but the SMC law produces redder spec- 
tra at shorter wavelengths for the same amount of extinction. 
Reddened templates with an SMC law reproduce quite well the 
optical spectra of dust- reddened quasars in the Sloan Digital 
Sky Survey (SP SS; seelRichards et alj|2003 ). while, using the 
Calzettil d 1997b law, IPolletta et all d2006h were able to repro- 
duce the SEDs of X-ray sources in the Spitzer SWIRE survey. 

The procedure of reddening a typical Type 1 quasar does 
not provide a satisfactory fit to the photometric data points of 
our sources. In Fig. [3] we show the results obtained when the 
prescription of iMignoli et al. (2004) for the extinction [SMC 
extinction law and E(B — V) = 0.7] is adopted. 

The dashed line shows the reddened quasar template nor- 
malized to fit the R-K s colour. For comparison, the unobscured 
quasar template is also shown (dotted line). Although the R-K s 
colour is obviously reproduced, the overall SED is not well re- 
produced, since the observed IRAC and 24 pm flux densities 
are systematically lower than predicted (up to a factor of 10 at 
24 pm). 

The discrepancy between the data points and the reddened 
Type 1 quasar template might be due to the application to ac- 
tive galactic nuclei of an extinction curve derived from galax- 
ies. The different behaviour of AGNs from galaxies can be 
attributed to different dust distribution (i.e., torus shape in 
AGNs), and gas-to-dust ratios, which can lead to an unusual 
dust reddening curve for AGNs. 



Since a reddened quasar template does not reproduce the 
shape of our data, w e adopted the red quasar template from 



Poll etta et al.l (120061) . which is a composite spectrum: in the 
optical/near-IR band, it is the spectrum of the red quasar 
FIRST J013435.7-093102 from|Gregj_etaD d2002|), while the 
average of several bright quasars from the Palomar-Green (PG) 
sample (Schmidt & Green 1983b with consistent optical data 
has been used in the IR. The IPolletta et al. ( 2006 ) template re- 
produces the observed data points significantly better, allowing 
for the observed sharp decrease from 0.2 to 0.7 nm i n the source 
rest frame (see Fig. [3] where the Polletta et alj|2006l spectrum is 
shown as a solid line). 

For source Abell 2690#29, which has a spectroscopic red- 
shift, the SED normalization has been obtained through a best- 
fit procedure. For source PKS 0312-77#45, where only a min- 
imum redshift was available prior to this analysis, we have left 
free to vary both the normalization and the redshift. 

In Table [3] we report the derived redshifts, column densi- 
ties Nh and de-absorbed L2-10 kev also for these sources which, 
similarly to the elliptical-like sources, belong to the Type 2 
quasar population. 

5. Physical parameters 

5.1. Bolometric correction 

Once the SED of the nuclear component has been determined, 
the following step is to estimate the bolometric luminosity Lboi, 
which is a quantity directly related to the central black hole 
activity. 

The bolometric luminosity L/, / can be estimated from 
the luminosity in a given band b, L/,, by applying a suit- 
able bolometric correction kboi,b - LboilLb- For the X-ray se- 
lected sources, the bolometric luminosity is typically estimated 
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from the luminosity in the 2-10 keV band feo/,2-10 kev - 
Lboi / L2-10 kev)- in previous work s, several authors used the 
bolometric correction obtained by Elvis et al.1 (1994) for lumi- 
nous, mostly nearby quasars, i.e., ^ o /,2-i0 kev~ 30. However, 
these corrections could be affected by the following uncer- 
tainties: firstly, they are average corrections obtained from 
a few dozens of b right quasars; secondly, as discussed in 
Marco ni et al. (2004), these corrections could overestimate the 
bolometric luminosities since they are based on the integral 
of the observed SEDs of bright unobscured AGN, without 
removing the IR bump (hence counting twice a fraction of 
~ 30% of the intrinsic optical-UV radiation). At lower lumi- 
nosities (typical of Seyfert galaxies, i.e., 10 42 - 10 44 erg s ), a 
lower value for this correction (kb i,2-io kev~ 10) was suggested 
(e.g.. iFabianl 120041) . For heavily obscured luminous sources, 
only few objects have been studied in detail; in particular, for 
two SWIRE Compt on-thick (i.e., log N H > 24 cirT 2 ) AGN 
Pollettaet al.1 (f2006) found kboifl3~& kev ~ 3 and « 100. 

In this work, thanks to the multi-band observations and ef- 
forts in disentangling the AGN and the host components, we 
try to derive directly the nuclear bolometric luminosity of our 
sources without assuming any average correction. We estimate 
Lboi by adding the X-ray luminosity integrated over the entire 
X-ray range (Lo.5-500 kev, not corrected for absorption) to the 
IR luminosity (Z4-1000 ,m)- 

io.5-500 kev has been estimated from the observed L2-10 kev 
luminosity assuming a single power-law spectru m with r=1.9 
(typic al for AGN emission; see, e.g., Fig. 6 of I Vignali et al 



tion has been performed in the 1-30 /mi range, after normal- 
izing the two SEDs to the same luminosity in the 30-100 yum 
range, where the anisotropy is thought to be negli gible. The 



2005 and references therein) plus absorption (where the col- 
umn densities are taken from the X-ray spectral analysis; see 
Perol a et al. I l2004l and Lanzuisi et al., in preparation) and an 
exponential cut-off at 200 keV. The median value found for the 
ratio Lo.5-500 kevlLz-io kev is ~ 4. The IR luminosity has been 
estimated by integrating the SED from 1 fim to 1000 /im using 
only the nuclear component for t he AGN hosted in t he ellipti- 
cal galaxies (see §4. Il l, and the iPolletta et alj (|2006) template 
for the point-like sources. 

Before computing the bolometric output of our sources, the 
derived IR luminosities must be properly corrected to account 
for the geometry of the torus and its orientation. The first cor- 
rection is related to the covering factor / (which represents the 
fraction of the primary optical-UV radiation intercepted by the 
torus), while the second correction is due to the anisotropy of 
the IR emission, wh ich is a function of the view ing angle (see 
Pier&Krolikll993l and lGranato & Danesell9"94l for further de- 
tails). 

We estimated the first correction (» 1.5) from the ratio 
of obscured (Compton thin + Compton thick) to unobscured 
quasars as require d by the most rec ent X-ray background syn- 
thesis model (see iGilli et al.ll2007l) in the luminosity range of 
our sources. A correction of « 1.5 implies an average covering 
factor which, in a simple torus geometry, corresponds 

to an angle 0*48° between the perpendicular to the equatorial 
plane and the edge of the torus. 

A first-order es timate of the anis otropy factor has been 
computed from the Silv a et al. (120041) templates as the ratio 
(R) of the luminosity of a face-on vs. an edge-on AGN, whose 
obscuration is parametrized as a function of Nh. The integra- 



derive d anisotropy factors are large only for the ISilva et al 
( 2004) template with higher column density (R ^3-4 for Nh — 
10 24 5 cm -2 ); since all of our targets are characterized by lower 
obscuration, such corrections do not affect our IR luminosities 
significantly (R * 1.2-1. 3 for N H ~ 10 22 - 10 23 4 cirT 2 ). In 
conclusion, the final combined corrections to be applied to the 
observed IR luminosities of our sources, given their column 
densities, are in the range w 1.8-2.0. After adding the X-ray lu- 
minosities, the IR correction factors would translate in a mean 
correction factor of ss 1.7 in the computation of the bolometric 
luminosities. 

In Table [3] the derived bolometric luminosities are reported 
along with the full range (Lo.5-500 kev) of X-ray luminosities, 
the IR (Li-1000 pan) luminosities and the bolometric corrections 
(hoi,2-w kev)- We note that our Li_iooo m estimates (hence L bo i) 
are robust despite the choic e of our SEDs. By comparing the 
£1-1000 nm obtained using the lSilva et al. (2004 ) model with the 
L\- 1 000 um obtained adop ting other recent average quasar SEDs 
(i.e. jRichards et al.ll2006l) . we have verified that the uncertain- 
ties in Li-iooo^im are within the a 10% level. The bolometric 
output of our targets is dominated by the IR reprocessed emis- 
sion, the primary X-ray radiation (Lo.5-500 kev) accounting only 
for <15% of the total luminosity. 

The derived median (mean) value of kb i,i-w kev is ~ 25 
(35+9; see Fig.|4]and colum n 8 of TableJUl, co nsistent with the 

1994 wid ely adopted 



30 from Elvis et 



value £fc /,2-io kev 
in pas t works. However, as pointed out also by lElvis et al 



(1 19941) . the bolometric corrections span a wide range of values 
(« 12-100); as a consequence, the adoption of a mean value 
could lead to inaccurate results. 

We note that in the extreme case where no corrections for 
the covering factor and the anisotropy of the torus are applied, 
we would obtain a median kboi,i-w kev of* 16. 

In Fig.[4]the derived Lb,,i as a function of L2-10 kev is shown; 
th e dot-dashed line join s the expected values from the analysis 



of iMarconi et al.l y004), where kb„i,2-\o kev is derived by con- 
structing an AGN refe rence template taking into account how 
the spectral index a ox (Zamorani et al. 119811) varies as a func- 
tion of the luminosity dVignali et alj|2003l) . Although the bolo- 
metric luminosities estimated for our o bjects are on average 
lower than those expected on the basis of lMarconi et al.l (|2004) 
relation, they are however consistent with a trend of higher 
khoi,2-io kev for objects with higher X-ray luminosity. If we fit 
our objects in the L\, i - -Li-iofey plane with the same slope as 
the IMarconi et al.l d2004l) relation, the difference in normaliza- 
tion is » 50%. 

5.2. Galaxy and black hole masses, and black hole 
Eddington ratios 

For the AGN hosted in elliptical galaxies we are able to infer 
both the galaxy and the blac k hole masses. The galaxy masses 
are estimated, assuming a ISalpeterl d 1955b initial mass func- 
tion (IMF), from the K s luminosities taking into account that 



Table 3. Inferred rest-frame properties of our targets 



Source Id. 


z" 


N b 

(10 22 cm- 2 ) 


^2-10 keV 

(10 44 ergs- 1 ) 


^0.5-500 keV 

(10 45 ergs' 1 ) 


-^1-1000 fan 

(10 45 ergs" 1 ) 


Lboi 
(10 45 ergs" 1 ) 


Lbol/Ll-lO keV 


4 

(10 n L,, o ) 


M star 
(10"M Q ) 


(10 9 M o ) 


(Lboi/ L E dd) h 


Abell 2690#75 


1 30 +,uo 


6.9 


3.2 


1.3 


9.7 


11.0 


34.7 


5.2 


3.5 


1.3 


0.065 


PKS 0312-77*36 


90+ 05 


1.0 


0.7 


0.2 


1.2 


1.5 


20.6 


1.0 


0.8 


0.2 


0.058 


PKS 0537-28*91 


1 30 +a4 ° 

• J -0.70 


25.8 


5.3 


2.6 


4.4 


7.1 


13.4 


2.8 


1.5 


0.7 


0.084 


PKS 0537-28*54 


>1.30 


1.6 


2.0 


0.6 


3.9 


4.6 


23.0 


3.0 


2.0 


0.7 


0.049 


PKS 0537-28*111 


1 20 +0 - 20 


9.1 


1.7 


0.7 


1.4 


2.1 


12.3 


7.8 


6.2 


2.1 


0.008 


Abell 2690*29 


2.08 


2.1 


8.4 


2.8 


78.4 


81.2 


97.0 


4.17 








PKS 0312-77*45 


1 8^+0-20 


8.0 


6.2 


2.6 


27.2 


29.8 


47.9 


1.65 








BPM 16274*69 


1.35 


2.5 


2.4 


0.8 


5.9 


6.7 


28.2 


8.8 


6.0 


2.5 


0.022 



" Photometric redshifts as derived from the analysis presente d in this paper (see 34.11 ). For source PKS 0537-28*54, only a minimun redshift was estimated; for sources Abell 2690*29 
and BPM 16274*69, the spectroscopic redshifts measured by Maiolino et al. (2006) are reported. 

b The column densities, measured through X-ray spectral fitting (see Perola et al. 2004 and Lanzuisi et al., in preparation, for details), were "matched" to the redshift used in the SED 
best-fitting procedure (using the relation N H (z)=N w (z = 0)x(l + z) 2 ' 6 ). 
c Absorption-corrected X-ray luminosity. 

d The 0.5-500 keV luminosities have been derived from the observed 2-10 keV luminosities as described in the text. The luminosities are not corrected for absorption. 

e The 1-10 00 luminosities have been derived from the integral of the nuclear SEDs (including the corrections described in the text). The values reported for the point-like sources 

refer to the lPolletta et alj |2006) red quasar template (see !4.2l for details). 

' The K,-band luminosities refer to the nuclear component for the point-like sources (Abell 2690*29 and PKS 0312-77*45) and to the host-galaxy starlight for the AGN hosted in the 
elliptical galaxies. 

g - h The reported M BH and Lboi/L E dd have been computed from the local L K - M BH relation jMarconi & Hunll2003h , under the hypothesis of an evolution of the M BH jM aar ratio of a factor 
two with redshift (see |5.2| for details). 
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Fig. 4. Bolometric luminosity vs. absorption-corrected 2- 
10 keV luminosity for the six AGN hosted in the elliptical 
galaxies (circles) and the two point-like AGN (triangles). The 
filled symbols refer to the values corrected for the covering fac- 
tor and torus anisotropy, while the empty symbols refer to the 
bolometric luminosity without applying these c orrections. The 
dot-da shed line represents the correlation from iMarconi et al 
(2004). The three dashed lines represent the loci of kf, o i )X =l0 
30 and 100. 
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Fig. 5. Black hole mass (Mbh) vs. bolometric luminosity (Ly) 
for the AGN hosted in the elliptical galaxies. The Eddington 
luminosity for a given Mbh is reported in the right-hand axis. 
The black hole ma sses have been estima ted from the local 
Lk - Mbh relation dMarconi & Huntll2003l) under two differ- 



ent hypotheses (see I5.21 >: (1) evolution by a factor of two of 
the MBH/M slar ratio with redshift in comparison to the local 
values (black filled circles); (2) no evolution of the -Mb#/M, M( . 
ratio (empty circles). The three dashed lines represent the loci 
of L bo i/L Edd = 0.01, 0.033 and 0.1 (from left to right). 



M star /LK for an old stellar population can vary from » 0.5 t o 
~ 0.9 (for ages between 3 and 6 Gvn lBruzual & Charlotl2003l) . 
We can derive M star directly from Lk since for these sources 
the ,-band emission is dominated by the galaxy starlight. The 



rest-frame Lk have been derived using the appropriate SED 
templates (see Sect. 14. \\ . The inferred stellar masses are in 
the range (0.8-6. 2)xlO n M , implying that our obscured AGN 
are hosted by massive elliptical galaxies at high redshifts. In 
Table |3] both the L K and M star v alues are report ed; we note 
that the different assumption of the Chabrie rl d2003l) IMF would 
produ ce a factor » 1,7 lower masses (Idi Serego Alighieri et al. 

To estimate the black hole masses, we take ad vantage of the 
local Mbh - Lk relation (IMarconi & Huntl l2003) which, taking 
into account the M star jhK values, is expression of the intrin- 
sic Mbh - M star relation. Given the challenging measurements 
of high-redshift black hole masses, the behaviour of this rela- 
tion with redshift is still matter of debate and different authors, 
using different techniques, have found different results. 



Woo et al] d2006h and IPeng et all d2006l) derive a signifi- 
cant evolution of the Mbh - M star relation with redshift, be- 
ing the MBH/M s , ar ratio larger, at high reds hift, up to a factor 



w 4 in comparison to the local value. In the I Woo et a l. (2006) 
analysis the discrepancy with re spect to the local value is al- 
ready present at z = 0.36, while IPeng et al. I d2006h find an av- 
erage Mbh I M star a factor >4 times larger than the local value 
at z > 1.7, while at lower redshifts (1 < z < 1.7) they de- 
rive a ratio which is at most two times higher than the local 
value, and may be consistent with m argin al or no evolution. On 
the other hand, IShields et aD d2006l) and iHopkins et all (§006) 
suggest that the MBH/M s , ar ratio is not significantly higher (at 
most a factor of two) than that measured locally up to z < 2. 

Given these uncertainties about the evolution of the Mbh - 
M s tar relation with redshift, we have estimated the black hole 
masses for our objects under two different hypotheses: (1) the 
MBH/M star ratio is higher than locally by a factor two in the 
redshift range (0.9 < z < 1 .4) of our sources; (2) the MBH/M slar 
ratio does not evolve with redshift. 

In both cases, our results imply very massive black hole 
masses (see Fig. O, in the range « 2.0 x 10 8 - 2.5 x 10 9 M in 
the former hyphothesis, with a factor of two lower values under 
the second hypothesis. 

Our estimated Mbh are consistent with the results derived 
bv lMcLure & Dunlopl (2004) studying a large sample of Type 1 
SDSS quasars and deriving the black hole masses from virial 
methods; most of their black hole masses are in the range 
1.5xl0 8 -2.5xl0 9 M in the redsh ift interval of our sample (see 
Fig. 1 of lMcLure & Dunlopll2004 . 

From the comparison of the bolometric luminosities com- 
puted in the previous section (see Table [3} with the Eddington 
luminosities calculated from the black hole masses esti- 
mated above, we derive that our obscured AGN are radiat- 
ing at a relatively low fraction of their Eddington luminos- 
ity (A » 0.008-0.084 and A sb 0.015-0.170 under the two hy- 
potheses; see Fig. and Table |3). This finding confirms and 



exten ds to a larger sample the results found by Maiolino et al. 
(2006) for two sources of our sample and by lBrusa et al.ld2005l) 



for a sample o f ERO s in the "Daddi Field". As suggested by 



Maiolino et al.1 d2006l) . the data indicate that our very massive 
black holes may have already passed their rapidly accreting 
phase and are reaching their final masses at low accretion rates. 
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Fig. 6. Bolometric luminosity as a fraction of the Eddington lu- 
minosity vs. redshi f t for th e whole sample of SDSS quasars of 
McLure & Dunlop (2004), plotted as small crosses. The large 
circles indicate the six HELLAS2XMM AGN hosted in ellipti- 
cal galaxies (symbols as in Fig . O ■ 



The estimated radiating efficiencies are sign ificantly lower 
than t he average Lb i/L E dd ~ 0-4 inferred by iMarconi et al 



(2004). However, since in the lMarconi et al.l d2004l) model only 



the phases of significant black hole growth are considered, our 
results are not in contrast with the proposed model but sug- 
gest that our targets belong to the tail of the sources charac- 
terized by low accretion rates. Consistently, our data (black 
filled and empty circles representing the evolution and no- 
evolution hypothesis, respectively) lie in th e lower envelope of 
the E ddington ratio distribution found by Mc Lure & Dunlopl 
(2004) br their large SDSS quasar sample. This is shown in 
Fig. [6] where our data are overlaid on the SDSS data points. 
This suggests that the SDSS quasar survey and the HELLAS 
survey probe different regimes of AGN activity: the SDSS sam- 
ples the brightest sources in the sky (R 1 20), most likely char- 
acterized by a high accretion rate, while our targets (X-ray se- 
lected, optically faint, i.e., R > 24, and obscured), are asso- 
ciated with a different evolutionary phase. We argue that the 
SMBH in our targets has already reached its final mass and the 
observed emission is witnessing a late stage of the accretion 
activity. 



6. Conclusions 

We have performed with Spitzer a pilot program to study a sam- 
ple of eight Type 2 (i.e., luminous and obscured) quasars at 
high redshift, selected from the HELLAS2XMM survey. Three 
sources have a measured spectroscopic redshift (two secure and 
one tentative) from near-IR spectroscopy; the remaining ob- 
jects have an estimated minimum redshift obtained from the 
R - K colours. On the basis of their /f s -band morphological 
properties, the sample is divided into two classes: sources with 
radial profiles typical of elliptical galaxies and point-like ob- 
jects. The most important results can be summarized as fol- 
lows: 



All of the eight sources have been clearly detected in both 
IRAC and MIPS 24 pm bands. 

The Spitzer observations have allowed us to detect the nu- 
clear component (often hidden at short wavelengths by the 
host galaxy) as thermal IR re-processed emission from the 
circumnuclear torus. While for the two point-like sources 
the nuclear component dominates at all frequencies, for the 
six sources with elliptical-like radial profile the contribu- 
tion from the strong stellar continuum is dominant up to 
the first IRAC bands, but the torus emission accounts for 
the entire emission at 24 pm. 

Taking advantage of the new Spitzer data, the nuclear SEDs 
of the sources have been modeled and new photometric 
redshifts have been estimated, following two approaches: 
for the elliptical sources, the nuclear emission has been 
"cleaned" from the host galaxy contribution adopting a 
two-component model (galaxy plus nuclear component), 
constrained using all the extensive observed data sets. For 
the point-like sources, the SEDs appear inconsistent with 
an extinguished Type 1 quasar template, bei ng well repro 



duced by an empirical SED of red quasars ( Pollett a et al 



2006). We find an overall agreement between the SED tem- 



plates and the data points, and the derived photometric red- 
shifts are consistent with the spectroscopic ones for two 
sources. 

Using the model components to extrapolate the nuclear 
SEDs in the far-IR regime, we derived the bolometric lumi- 
nosities (being in the range « 10 45 - 10 47 erg s -1 ) by adding 
the IR luminosities to the full range of X-ray luminosities. 
In this computation, we have considered and discussed the 
corrections to be applied to the observed IR luminosities to 
take into account the covering factor of the torus and the 
anisotropy of the IR emission. The median 2-10 keV bolo- 
metric correction is « 25, consistent with the value typically 
assumed in literature. 

For the elliptical sources, thanks to the independent esti- 
mates of the stellar light and nuclear bolometric luminos- 
ity, the physical parameters of the central black holes have 
been estimated using the M BH - L K relation and exploring 
different hypotheses for the evolution of the M BH /M star ra- 
tio with redshift. Under the hyphothesis that the MBH/M slar 
ratio is a factor of two higher at z ~ 1.2 than locally, our 
luminous, obscured AGN have masses in the range (0.2- 
2.5)xl0 9 M , reside in massive [(0.8-6.2)xlO u M ] high- 
redshift ellipticals and are characterized by low Eddington 
ratios (Ax 0.008-0.084). Through our direct estimate of the 



IR luminosity, we confirm the conclusion of Mai olino et al 



(2006) that these black holes may have already passed their 
rapidly accretion phase. 
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